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ing wound string interactions in an impact parameter picture, we show that a strong 
exponential suppression in the interaction rates for d > 3 spatial dimensions reflects 
the classical argument that string worldsheets generically intersect in at most four 
spacetime dimensions. This description is appropriate in the early universe if wound 
strings are heavy - wrapping long cycles - and diluted. We consider the dynamics of 
a string gas coupled to dilaton-gravity and find that a) for any number of dimensions 
the universe generically stays trapped in the Hagedorn regime and b) if the universe 
fluctuates to a radiation regime any residual winding modes are diluted enough so 
that they freeze-out in > 3 large dimensions while they generically annihilate for 
d = 3. In this sense the Brandenberger-Vafa mechanism is operative. 
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1 Introduction 



One of the few mechanisms aiming to explain the hierarchy between three large and 
six small spatial dimensions within superstring theory is due to a suggestion, some two 
decades ago, by Brandenberger and Vafa [Tj; see [2] for a review. In this scenario the 
early universe consists of a hot string gas in thermal equilibrium near the Hagedorn 
temperature. The topology of space has non-trivial cycles supporting winding modes 
in the gas. The background metric and string coupling evolve with the low energy 
effective dilaton-gravity equations of motion according to which the winding modes 
resist the expansion of the spatial directions they wrap. If due to a thermal fluctuation 
a number of dimensions starts growing then eventually the equilibrium number of 
winding modes will drop to zero. The winding modes have the capacity to relax 
to equilibrium through annihilations with anti-winding modes; if these interactions 
are efficient then at large volumes the winding numbers will vanish allowing the 
corresponding dimensions to grow. The Brandenberger- Vafa (BV) mechanism relied 
on a simple dimension-counting argument that wound strings generically intersect in 
at most 3 spatial dimensions, singling this out as the maximum number of dimensions 
in which winding numbers have the capacity to track their equilibrium values, thereby 
dropping to zero and allowing the dimensions to grow large. In [3] this argument 
was supported using numerical simulations of a network of classical strings, though 
gravitational dynamics was not taken into account. 

Over time it became clear that strings are not the only fundamental degrees of 
freedom of string theory and that higher dimensional objects (membranes) are also 
fundamental states of the theory; superstring theory was shown to result from com- 
pactification of a higher dimensional theory, M-theory. In a paper by Alexander, 
Brandenberger and Easson [Ij the setup of [I] was extended to include p-branes for 
p = 0,1, 2, 4, 5, 6, 8 in the weak-coupling limit of M-theory with one small dimension 
compactified on 5*^ (type IIA string theory). The other spatial dimensions were com- 
pactified on a 9-torus. The authors argued that fundamental string winding modes 
are still the decisive objects regarding decompactification and that the conclusions 
of still hold. They also pointed out that a further hierarchy between dimensions 
could arise. Past the string scale, as the universe grows, more and more energy is 
needed to support wound branes of highest p, hence highest p branes would tend to 
decay first. As two p-branes can intersect in at most 2p+l spatial dimensions, there 
is no obstacle for the disappearance of p-branes for p > 2. But 2-branes can allow for 
a 5-dimensional subspace to grow first. Further, within this subspace, 1-branes will 
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only allow for a 3- dimensional space to continue expanding, as in [T] , hence one is left 
with a 3-2-4 dimensional hierarchy. 

These claims relied on heuristic thermodynamic and topological arguments. Aim- 
ing to carry out a more rigorous investigation, Easther et al. ^ considered the full 
equations of motion for IID supergravity on a homogeneous but anisotropic toroidal 
background, coupled to a gas of branes and supergravity particles. Focusing on the 
late time behavior of the system, they justifiably ignored excitations on the branes 
and included only M2-branes, since M5-branes (the other fundamental states of M- 
theory) would annihilate efficiently in the full 11- dimensional spacetime. Motivated 
by the BV mechanism and the arguments of the authors of [5] chose initial states 
resulting from fluctuations that would leave 3 dimensions unwrapped, some num- 
ber of dimensions partially wrapped and some fully wrapped. The conclusion was 
that indeed the dynamics leverage the topological reasoning and a hierarchy among 
dimensions is established. This conclusion was supported further by [B], where in 
addition nontrivial fluxes were included (the 3-form gauge fleld of IID supergravity). 
In fact, the presence of fluxes seemed to enlarge the possible space of initial condi- 
tions that lead to three large dimensions at late times. Speciflcally, for the case of six 
initially unwrapped dimensions, the dynamics of fluxes introduced a new hierarchy 
suppressing the growth of 3 out of the 6 unwrapped dimensions. 

An apparent limitation of the BV argument is that it seems to depend crucially 
on non-contractible spatial cycles and their associated topologically stable winding 
modes. Phenomenologically viable compactifications of string theory, however, may 
not have such cycles. Nonetheless the authors of [7] surmised that these more general 
spaces might still support "pseudo-wound" modes, long strings that extend around a 
dimension but are contractible. If these strings are stable over time scales larger than 
the cosmological Hubble scale, then as far as the dynamics are concerned they play 
the same role as stable wound strings. In [7], using numerical simulations for string 
networks on toroidal orbifolds with trivial fundamental group, the authors showed 
that pseudo-wound strings generically do persist for many Hubble times, suggesting 
that the requirement of non-contractible cycles can be relaxed. 

The results up to this point seemed promising, but it remained to actually test 
the heart of the argument: whether at early times, thermal fluctuations near the 
Hagedorn era and string (or brane) interactions really lead to annihilation of winding 
modes in a 3-dimensional subspace. An early attempt to investigate this was carried 
out in [8J. The authors considered a gas of 2-branes and supergravity particles, along 



2 



with excitations on the branes that lead to a hmiting Hagedorn temperature. This 
setup was within the low-energy limit of M-theory compactified on a 10-torus, with an 
anisotropic and homogeneous metric evolving according to 11-dimensional Einstein 
gravity. The winding numbers of 2-branes evolved according to Boltzmann equations. 
The authors assumed initial conditions in which the total volume of the torus was 
fixed but otherwise assumed that all states were equally likely. By numerically solving 
the coupled Boltzmann-gravity equations the authors concluded that the number of 
unwrapped dimensions at late times depended crucially on the initial volume of the 
torus. Typically a large (and monotonically increasing) overall volume would decrease 
the interaction cross-section of branes too quickly, eventually leading to brane number 
freeze-out. If the initial volume was constrained according to holographic arguments, 
the initial winding numbers proved so small that all dimensions would decompactify 
early on. Three dimensions was not found to be singled out by the dynamics. 

Similar all-or-nothing behaviour was found in [21 [ID] for IIA theory compactified 
on T^. In these papers the dilaton-gravity equations for the background were coupled 
to Boltzmann equations for winding modes and radiation. Even though this behaviour 
was attributed to the rolling of the dimensionally reduced coupling to weaker values, 
we emphasize here a more decisive phenomenon that prevents the annihilation of the 
winding modes and yields the all-or-nothing behaviour. If the initial energy density 
of the universe is large, the system is found in the Hagedorn phase with a significant 
amount of winding present in thermal abundance but in a regime that resembles a 
matter-dominated universe with vanishing pressure. This, along with "friction" due 
to the dilaton's velocity, results in an insignificant growth of the wrapped dimensions 
(even over an infinite amount of time). With the total energy nearly constant, the 
equilibrium number of winding modes generically does not fall to zero. In this sense 
the system stays "trapped" in the Hagedorn phase. It is very likely, however, that this 
problem is particular to the approximation of treating the background with the lowest 
order dilaton-gravity dynamics. Corrections to these ([20]), or a different treatment 
of the metric degrees of freedom, could alleviate it. An alternative approach would 
be to keep the background dynamics to lowest order and still consider a high density 
initial phase - a fairly natural assumption - but consider large volume fluctuations 
that could yield an exit. This is the approach we adopt here. 

Finally, and most importantly for our purposes, the aforementioned problems 
were independent of the number of dimensions growing large. The reason was that 
the rate at which wound strings annihilated only fell off like the inverse volume of 
the transverse dimensions. This failed to single out three large dimensions as special. 
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suggesting that the Brandenberger-Vafa argument might not be supported by the 
dynamics underlying string/M-theory. 

In this paper, we re-examine this conclusion and suggest a possible way in which 
string dynamics may indeed favor three large dimensions. Our basic approach is this: 
According to the Brandenberger-Vafa dimension-counting argument, one expects that 
string interaction rates should be dramatically suppressed when the number of large 
spatial dimensions is bigger than three. Moreover, as the dimension-counting argu- 
ment is purely classical, one expects it to be valid in a regime where the wound strings 
behave nearly classically and can be regarded as one-dimensional extended objects 
tracing a two-dimensional worldvolume. In such a regime the quantum thickness of 
the strings should be small compared to their length along the dimension they wrap 
and also small compared to the size of the transverse space. This suggests, in contrast 
to our previous work, that we hope for a dilute gas of winding strings. Furthermore, 
as we discuss, in a dilute regime we are lead to work in an impact parameter repre- 
sentation of the string scattering amplitude. As we will see, this makes manifest the 
distinction between three and more large spatial dimensions regarding the interac- 
tions of winding modes. Our main observation is that if the universe fluctuates out of 
an initial dense Hagedorn regime - something that we believe is generically necessary 
in order to match to a realistic expanding cosmology in the string gas scenario inde- 
pendently of the decompactification mechanism - then any residual winding modes 
that were thermally excited in the Hagedorn phase, are indeed diluted enough so that 
they freeze-out in c? > 3 spatial dimensions. Further, for d = 3, the enhancement in 
interaction rates due to the length of wound strings generically overcomes the suppres- 
sion due to the weak coupling and winding modes may annihilate efficiently. While 
a number of important issues remain, this appears to be the first demonstration of 
dynamical string theory decompactification that generically yields three large spatial 
dimensions. 

By way of outline we begin with a discussion of the impact parameter represen- 
tation, proceed to set up our model for the string gas, and finish with a numerical 
simulation, along the lines of [9], that will allow us to identify the regions of phase 
space in which three or more spatial dimensions decompactify. 
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2 Interaction amplitudes and impact parameter 
picture 



In this section we derive interaction rates for wound strings in a semiclassical impact 
parameter picture. We will show that when long strings interact at impact parameters 
larger than their thickness, there is an exponential suppression in the interaction rates 
for d > 3. 

The starting point is the Virasoro-Shapiro amplitude for wound strings in ci = 
D — 1 large dimensions given by [HI [T2] 



with s computed either from the right-moving or left-moving momenta of the closed 
string, s ~ AE? /a'"^ with R the radius of the dimension that the strings wrap. The 
imaginary part of the amplitude as t — i- is 



Here k|)_2 = n"^ /V is the gravitational coupling in D — 2 dimensions, where V is the 
transverse compactification volume times the area of the torus wrapped by the strings 
pT] . By the optical theorem ^Im(A(s, )f:))|t=o ~ a' k\,_2S controls string interactions. 
It is crucial to observe that for D = A this quantity is dimensionless and gives the 
probability for two colliding winding strings to interconnect and unwind (to leading 
order), while for D > 4 it has units of (length)^"'^ and represents a cross section in 
the D — A dimensions transverse to the moving strings. This reflects the fact that 
long strings generically intersect in D = 4, like point particles moving on a line, while 
they generically miss in D > 4 and the relevant quantity becomes a cross-section. 

One can consider the interaction probability in an impact parameter picture. As 
discussed above, long wound strings have an effective impact parameter in the D — 4 
directions transverse to the motion of both strings. The impact parameter h is the 
conjugate variable to the transverse momentum q = \f—i and the amplitude in this 
representation is obtained by the following transform in the transverse directions. 





ImA(s,t = 0) = — 4_2s' 



(2) 




(3) 
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For the Virasoro- Shapiro amphtude ([l]), using g ^ = «' dxx°'''^^ ^, this gives 

where Y = log(^) and 7(0, x) is the lower incomplete gamma function. The imagi- 
nary part of the above amplitude in the limit ^ Ya' is 

T At h\ . 7^a'Kl_^S 

Imylls, 0) — )■ — — — — .^.^ 4^0' (5) 



These results are similar to those found in fT3], the difference being that the au- 
thors of [13] consider graviton scattering and take the number of transverse directions 
to be D — 2. In fact the interpretation of 6 as a classical impact parameter in ([s]) can 
be justified along the lines of [13] • In the high energy limit (s — )■ 00 which for wound 
strings is — ?■ 00 - precisely our limit of interest) the authors of [I3J sum up the am- 
plitude to all loop orders to a unitary eikonal form. The large R or large energy limit 
localizes strings in the transverse directions and reveals classical behaviour, much as 
the eikonal treatment in quantum mechanics (or optics) reveals semiclassical particle 
(or ray) behavior. 

Note that A(s, h) is dimensionless for any D. It determines the annihilation prob- 
ability P{h) via 

Pih) = -\miAisM (6) 

V 

with lmA{s^ h) as in equation ([s]) and v the velocity of the colliding strings in their 
center of mass frame. This prescription can be shown to satisfy the usual unitarity 
conditions in the large s limit [MJ 115] . 

The quantity Ax^ = AY a' = 4a' log (i?^/ a') appearing in ([s]) is interpreted as 
the quantum thickness of the string. It measures the fluctuations about the classi- 
cal straight string configuration. The fact that it increases logarithmically with the 
string's length reflects that it is energetically less costly to excite oscillators on a long 
string. Similar string spreading effects occur in high energy collisions and for strings 
falling into black holes [TU]. Note that this string spreading does not include the 
effect of real (as opposed to virtual) oscillator excitations as would be appropriate in 
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the Hagedorn phase of a string gas. In the Hagedorn phase wound strings are highly 
excited and their spread in the transverse directions is comparable to the length of 
the dimension they wrap [IZ!]^ These wiggly strings are very likely to intersect, lead- 
ing to rapid interactions which keep the strings in equilibrium. But as the universe 
expands and cools down the equilibrium phase becomes one of pure radiation. Then 
the oscillator excitations decay away and the spread of the wound strings approaches 
Ax. This justifies our use of the amplitude (|5]) if 6 > Ax. 

It is useful to contrast the impact parameter picture to the more standard method 
of obtaining a scattering probability. Typically one derives a cross section a and the 
collision probability is simply na where n is the number of targets per transverse 
volume. If one has a collision probability in impact parameter space, P{h), then the 
scattering cross section is obtained via [HI [19] 



In other words, na is an averaged probability in impact parameter space. Most 
frequently it is assumed that the targets are uniformly distributed {n = constant) 
and one obtains 



In the case of the optical theorem, for example, one can immediately derive a = 
^lmA{s,t = 0) using equations (jsj), ^ and ([3|. It is the n = constant assumption 
that we are willing to relax here. There are two ways in which it can be justified. 
First, if the targets are dense and uniform as in collider experiments. A test particle 
in that case will interact with targets at all impact parameters so one can integrate as 
in Q. Second, if the time between collisions is much smaller than the total time over 
which collisions take place. Then the test particle is given enough time to interact 
with targets at all impact parameters (assuming each collision is at a random impact 
parameter) and the averaging over impact parameters is essentially a time average. 

But if the winding modes in a string gas are dilute, with a mean separation much 
larger than their thickness, the dense target assumption above does not apply. It 
could still be that, since the strings move in a compact space, they collide repeatedly 
with each other and a time average is appropriate. It then becomes a matter of 

^In the Hagedorn phase strings perform a random walk in aU directions. As their energy scales 
with their length, their mean extent in all directions scales as \/E. This is the dependence of the 
winding number on energy as we will see in the section on thermodynamics. 




(7) 
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timescales. We need to compare the mean time between collisions with the recoUapse 
time, the time required for winding modes to pull the universe back to a small-radius 
regime where winding modes are no longer dilute. An additional effect which must 
be taken into account is that the string coupling is time dependent. This could also 
invalidate the use of a time-averaged cross section. 

We thus have to develop a model for the distribution of interactions over impact 
parameters. We will return to this in the next section after we set up the rest of the 
dynamics. 



3 Equations of motion 

In this section we write down coupled dilaton-gravity and Boltzmann equations for 
the matter degrees of freedom. For further details on the thermodynamic phases and 
energy conservation see [9l [20] and references therein. 

The general setup is similar to [9], except that instead of an anisotropic universe 
we consider d growing dimensions all with the same radius, and hold the remaining 
9 — d dimensions frozen at the self-dual radius. By removing the randomness in the 
choice of initial radii present in [0] we can see more clearly the dependence of the 
winding annihilations on the number of growing dimensions. 

We consider type IIA string theory with a flat FRW metric on a torus for the 
d = D — 1 growing dimensions, 

ds^ = -dt^ + a'e^^^^^J2idx^i < a;^ < 27r (9) 

and a homogeneous shifted dilaton <f{t). From now on we set a' = 1. When the 
metric and dilaton are coupled to matter, the equations of motion are 

^ = h^' + dx') 

2 (10) 
A = + —^e^P 

and the Hamiltonian constraint (Friedmann equation) is 

E = (27r) V^(0' - rfA') (11) 

Here E is the total energy in the string gas and P is the pressure (times the volume) 
in d dimensions. 
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3.1 Matter content and Boltzmann equations 

The background equations of motion are coupled to phenomenological Boltzmann 
equations that govern the evolution of matter. We model matter with three species: 

• Winding modes that evolve according to 

W = -TwiW^ - (Wf) (12) 

We specify the interaction rates F and the thermal equilibrium values (■) below. 
The total energy in winding and anti-winding modes is E]y = 2dWe^ and their 
contribution to the pressure is Pw = —2We^. 

• Radiation, or pure Kaluza-Klein modes, evolving according to 

k = -Tk{K' - {Kf) (13) 

The energy in KK and anti-KK modes is = 2dKe~^ and their pressure is Pk = 
2Ke-^. 

• Finally we include string oscillators, or massive string modes, as pressureless 
matter. The oscillator modes fill up the energy budget via Eo^c = E — {Eyy + E^). We 
do not need a Boltzmann equation for these modes since the dilaton-gravity equations 
of motion automatically conserve energy, dE = —PdV. 

3.2 Thermodynamic phases and interaction rates 

Near the self-dual radius the gas of strings is in a high density Hagedorn phase. The 
thermodynamics of this phase has been studied in [211 122] . The quantities of interest 
here are the equilibrium values of the winding and KK number^ 




Since E ^ 1 most of the energy in the Hagedorn phase resides in oscillator modes 

(£;osc^^- v^). 

^These values are derived under the assumption that the microcanonical energy is spht equally 
amongst all dimensions, see [S]- 
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As the volume of the universe grows and the energy density drops, the equihbrium 
state should be one with only radiation]^ As the condition for that transition we set 

^ < c,Ti^' (15) 

with Th the Hagedorn temperature, Vd = (27r)'^e'^'^ the volume of the large dimensions, 
and Cd a Stefan-Boltzmann constant appropriate to the IIA gas of 128 massless Bose 
and Fermi degrees of freedom in d dimensions. 

In this phase the equilibrium values are 

{W) = {K) = ^Ee' (17) 

That is, at equilibrium all the energy is in radiation (KK and anti-KK modes). 

Now we need to specify the interaction rates entering the Boltzmann equations. 
Recall that for winding modes, with an impact parameter b in D — 4 dimensions, 
the interaction probability is P{b) = ^lmA{s,b) with lmA{s,b) given in ([5|. For 
two wound strings moving in the Xi direction and with opposite winding along X2, 
the right moving momenta are p_Ri,2 = {E,±Ev,±R/a') so sr = —{pm + Pr2)^ = 
AE'^ ^ {2R/a'Y for slowly moving strings. Putting things together, the interaction 
probability per unit time (per winding mode in the direction of motion) can be written 
as 



w 



no 

4 V \a' I 



X 




with V the total spatial volume of the 9-torus. In terms of our variables, and with 

2 

a' = 1, we have ^ = 2^^"^ ^^"^ R = e^. Note that Fq is the interaction rate used 
in [9]. 

As explained earlier, an impact parameter representation is only appropriate in 
the radiation phase, when the separation between winding modes r is larger than the 



■^See [211 El] for a detailed treatment of the conditions for equilibrium between massive and 
massless modes in a string gas. 
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string thickness Ax. From the thermodynamic distributions of [22] we can estimate 
how the mean velocity w of a single winding mode depends on R and E (see appendix 
A)gThe mean time between collisions, or re-collision time, is then tr — ^- In practice, 
as we numerically integrate the equations of motion, once we are in the dilute regime 
we randomly choose an impact parameter b on every re-collision time. The impact 
parameter is chosen at random, from a uniform distribution in the transverse D — A 
dimensions, up to the maximum value b = r. In effect, under the assumption of 
isotropy and even distribution of winding modes, we treat their interactions as d- 
many copies of a single interaction over a periodic lattice r'^ {d = D — 1). We have 
now enough information to determine Th. 

Another concern, raised earlier, is that in the dilute regime the winding strings 
might not have time to collide before the universe re-collapses to a dense Hagedorn 
phase. In principle this could happen even in D = 4. We test for this as follows. Upon 
entering the dilute regime we estimate the re-collision time and turn off interactions, 
i.e. set r = 0. If the negative pressure from the frozen winding modes recollapses 
the universe in a time smaller than tr it means that freezing the interactions was 
consistent, that is, the winding modes truly had no time to collide. On the other 
hand, if after time tj. we are still in the dilute regime, then string interactions must 
be taken into account. 



Thus (18) provides our description of string interactions in the dilute regime. In 
the Hagedorn phase the strings have highly excited oscillator modes which enhance 
the interaction rates since more string is available. This was studied in [9] and it 
amounts to inserting an overall factor of yi? in the Boltzmann equations. 

We also need to specify the interaction rate for KK modes. Since the wavelength 
of these modes grows with R a semiclassical impact parameter picture at large R is 
not appropriate. Instead we should average over impact parameters. Since we already 
know the averaged interaction rate Fq for winding modes, by T-duality we can take 

= ro|A-s>-A- 



^Even though we are working off-equihbrium we consider the equihbrium velocities to be a good 
approximation. In other words, we are assuming kinetic equihbrium and explore the possibility of 
chemical equilibrium. To be precise, note that v — Wmis — \J (w^) is the root-mean squared velocity. 



11 



3.3 Initial conditions 



We need to integrate the coupled equations (10), (12), (13) subject to the constraint 
( |ll| . We need 6 initial conditions: 9?o, 00, Aq, Aq, Wq and Kq. In this section we 
describe our method for sampling from the space of initial conditions. The basic idea 
is to fix the value of 0, scan over allowed values of Aq and (po, and average over the 
values of Aq, Wq and Kq using a suitable probability distribution. 

For effective supergravity to be valid and to ensure weak coupling we take ipo = —1 
and ipo < —1. Recall that the dilaton evolves monotonically to weaker coupling with 
the absolute value of its velocity decreasing p. 

With 00 and ipo fixed we consider the universe at equilibrium with A = and 
A = (at the self dual radius). Calculating the energy E^q = E\j^^Q determines the 
equilibrium number of winding and KK modes per dimension the universe would have 
at the self dual radius, 

12 V TT ■ 

We will use these values to set upper (lower) bounds for choosing Wq (Kq) below, 
once the volume fiuctuates to larger values. 



Next we choose the scale factor Aq of the d-torusj^ We take Aq > since by 
T-duality we need not consider smaller volumes. Having already determined Eg^, 
choosing Aq fixes the energy density and thus the equilibrium thermodynamic phase. 
The choice of Aq can be thought of as merely a choice of initial condition but since we 
choose the winding and KK modes with respect to their values at A = (see below) 
it is more appropriate to interpret it as a fiuctuation from A = to a larger volume. 
Since we are choosing Aq at random we are then assuming the thermal distribution to 
be fiat. This is not an arbitrary choice. The entropy to leading order in the Hagedorn 
phase is 5*0 — E/Th and thus it costs no entropy or energy for such a fiuctuation. To 
next order the dependence of the entropy on the radii R = for d spatial dimensions 

is [221 [22] 

Si ~ log[l - T{2dY\'nEf'^-^e-'^^] 

with 7] ~ 1/ E?. Even though this contribution is very small for the ranges of energies 
and radii we consider, it would be interesting to study more precisely the effect that 
these corrections give to our scenario. 



^We consider the dynamics of d dimensions, while the remaining 9 — d are kept frozen at the self 
dual radius, A = with A = for all times. 
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The Hubble rate A can also fluctuate away from zero, which is the value that 
maximizes the entropy. In the Hagedorn phase the entropy is given to a good approx- 
imation by = E/Th- Thus Aq is chosen randomly from the Gaussian distribution 

oc e-^o/(2-|,) (19) 
with ajj = 2{2w)°d ■ radiation phase the entropy is = ^qV^ ( ^r^j 



Using (11), to leading order in Aq we have the distribution 

oc e-^o/(2-?) (20) 



with 



1 



a'r={^y''< P=^^ PH^c,Ttt' (21) 



It remains to choose the initial winding and KK numbers. Depending on whether 
we are in the radiation or Hagedorn phase, the equilibrium number of winding modes 
could be zero or not. Since we do not want to begin with zero winding (we would 
not be testing the BV mechanism in that case) the lowest value of W we may pick 
is 0.5, our chosen threshold between zero and non-zero winding. The furthest we can 
fluctuate from equilibrium (the largest winding) is W^d- However it's possible that 
the volume is so large that there isn't enough energy to support that much winding. 
This occurs if Wsd > Ee~^/{2d). Putting everything together, the initial winding 
number is chosen randomly in the range 

(^Max{0.5,(l^)},Min{(l^),d,^e-^}^ (22) 

In the Hagedorn phase the KK number can fluctuate between {K)sd and the equi- 
librium value at the given Aq, so we choose a value randomly in this range. In the 
radiation phase, given Wq, we compute the energy in winding Ew = 2dWoe'^. The 
rest of the energy should be available to radiation with the maximum KK number 
being -ft'max = {E — Ey/)e^ / {2d) . Therefore Kq is chosen randomly in the range 

((ir)sd, ^max)- If ^max < (^)sd We SCt Kq = K„ 



^When integrating the equations of motion we need to be careful not to produce more radiation 
than energy conservation allows. If at some point E = + Ey/, that is all the oscillators decay, 
we set {K) ^ {E - Ew)e^/{2d). 
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Once initial conditions are fixed we integrate the equations of motion until either 
the winding modes annihilate {W < 0.5) or the interactions freeze out, which we 
define as TqW < O.IH. We use the maximum rate Tq instead of the total Tw to allow 
for the possibility that, depending on the randomly chosen value for b, strings could 
interact even for D > 4. 



4 Results 



With 00 fixed at -1 we can scan initial conditions over a two-dimensional lattice of 
points (Ao, v^o)- For each lattice point we do 1000 runs to average over different choices 
of Ao, Wo and Kq. 

The results for values of d ranging from 9 to 3 are contrasted in figures [T] and |2} For 
all values of d we find that equilibrium is maintained during the Hagedorn phase. A 
consequence of this is that for a large range of initial conditions, as long as the system 
starts out in the Hagedorn phase, it will remain forever trapped there (see figure [2]). 
This is related to the limiting value of A(t) as t — )■ oo in the solution to the equations 



of motion (10) in thermal equilibrium [20] and it is also the same behaviour found in 
[23] at exact equilibrium. When the volume is large enough such that (W) (Qyet 
the system is still in the Hagedorn phase, then the universe decompactifies for any 
d. This region gets narrower as d increases as seen in figure [2j But for d > 3 this 
is essentially the only region in parameter space where decompactification occurs. 
If the system gets to the radiation the oscillators decay to massless modes 

and the long strings start diluting, hardly any choice of initial conditions leads to 
decompactification. For d > 3 and outside the Hagedorn phase there are few cases 
(of order 1%) in which the universe decompactifies. Those rare cases have very small 
winding number and happen to have collisions at small impact parameters. 

By contrast, for the d = 3 case shown in figure [T| we see that even in the radiation 
phase long strings are able to annihilate. It is interesting to note that the effects of 
large A enter in two competing ways. First, due to the factor of s ~ i?^ = e^^ in 
the amplitude, long strings interact more efficiently, even at weak coupling. But at 
large A the effect of dilution is more dramatic and strongly (exponentially) suppresses 
interactions for > 3. 

^In practice this would be (W) < 0.5. We round W to zero and not (W). 
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A comment on the choice 00 = —1. We could have considered smaller (absolute) 
values for ipo, still valid within supergravity. This can be compensated by a small 
(logarithmic) shift in the initial dilaton to smaller values such that the initial energy 
remains the same. The qualitative results should be unaltered. 



5 Summary and discussion 

The Brandenberger-Vafa mechanism relies on a classical dimension- counting argu- 
ment, namely that the worldvolumes of one-dimensional objects will generically in- 
tersect in at most three spatial dimensions. To see the mechanism at work one needs 
to be in a regime where strings behave as semiclassical one-dimensional objects. That 
is, one needs their length to be much larger than their effective quantum thickness, 
and also their thickness to be much smaller than the size of the transverse directions 
(the space is compact). We realized these conditions in a simple isotropic setup where 
the length of the string R was the same as the size of the compactification manifold. 
With oscillators excited, as in the Hagedorn phase, strings have a significant spread 
in all directions and the classical picture fails. But when oscillators decay, as in the 
radiation phase, the thickness of strings grows as Ax ~ y/log R. Thus strings begin 
to behave classically as R/ y/\ogR grows. To model string interactions in this regime 
we developed an impact parameter representation of the string scattering amplitude. 
This allowed us to show that in this regime the BV mechanism indeed operates and 
favors decompactification of three spatial dimensions. 

To enter this regime we had to consider departures from equilibrium, often large. 
Clearly in the radiation phase this is necessary since the equilibrium number of wind- 
ing strings is zero. In the Hagedorn phase strings rapidly come to equilibrium and 
the pressure vanishes. This means the universe tends to remain stuck in the Hage- 
dorn phase, and for some number of dimensions to decompactify a large fluctuation 
is needed, either in the Hubble rate or in the initial volume, to send the system to 



a regime where the equilibrium winding number is zero. As the distribution (19) 
typically allows for only small fluctuations in the Hubble rate, we had to consider 
large fluctuations in the volume to realize the BV mechanism. An important next 
step would be to understand the likelihood of such a fluctuation taking place in the 
early universe. 

One shortcoming of our framework was that, even though we considered a dilute 
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d=9 d=8 d=7 




d=6 d=5 d=4 




Figure 1: Number of cases decompactifying as a function of (po and Aq for different 
choices of growing dimensions d. For d = 4, 9 the 2- axis is chpped at 200 to make 
the fewer decompactifying cases visible. 
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1 2 3 4 5 



Ao 

Figure 2: A plot of the (Aq, v'o) plane contrasting the cases d = 3, 4, 6, 9. If the initial 
equilibrium winding number in the Hagedorn phase is non-zero the system typically 
stays trapped in the Hagedorn phase, unless the initial Hubble rate is large and the 
initial winding number is small (thin dark region labeled 'Equilibrate out of Hagedorn 
for Ao = lOo"'). If the initial equilibrium winding number is zero in the Hagedorn phase 
then the universe typically decompactifies in any number of dimensions (regions in 
the upper left corner, labeled by dimension). But if the universe begins in a radiation 
phase with a dilute gas of winding strings then only d = 3 will decompactify (orange 
region to the right of the grey line). 
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gas of winding strings, we modeled the resulting pressure as a homogeneous term 
in the gravity equations of motion. This led us to consider their back-reaction on 
spacetime in an all-or-nothing manner, in which any amount of winding would op- 
pose expansion while zero winding would not. As far as testing the interactions and 
eventual annihilation of strings, which was our focus, this shouldn't be a concern. 
But a more detailed investigation of string gas cosmology should address the issue of 
spatial inhomogeneity. Finally it would be interesting to extend the analysis of this 
paper to the more general context of M-theory, taking into account the effects of the 
full p-brane spectrum. 
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A Root-mean-square velocity of winding modes 

Thermodynamic quantities can be calculated using string distributions derived by 
Deo, Jain and Tan [211 l22] in the microcanonical ensemble. They show that the 
average number of strings with winding charge vector w, Kaluza-Klein charge vector 
k and energy e on a D-torus with total energy E is given by 

D{e, W, k, E) = ^Mi5g-«q^A-lq/4 ^23) 

e 

where 

E 



u = 


e{E - e) 


q = 


(w,k) 


N = 




\/detA 


A = 













One can consider a unit winding mode, wi = (1,0,...) along one of the d large 
dimensions, {Ri = R, i = 1, d and Rj = 1, j = d + 1, -D) and calculate the mean 
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momentum squared of the string. 



2^ ^ /o''&/ci^k e D(e,w = Wi,k,E) 
Jf^dejd^k D{e,w = wi,k,E) 

(24) 

~ 27r2 de^D/2Q-UTTR'^ 
Jo e 

Given the total energy E and the radius R the two integrals above can be evaluated 
with saddle point methods or numerically. For a heavy winding mode (large R) we 
have (f ^) ~ ^ and Wrms = a/T^- 
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